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Abstract 



Aims. Supersonic turbulence is a large reservoir of suprathermal energy in the interstellar medium. Its dissipation, because it is 
intermittent in space and time, can deeply modify the chemistry of the gas. This is clearly seen in the framework of shock chemistry. 
Intense turbulent dissipation also occurs in regions of large velocity shears, sharing with shocks the property of intermittency. Whether 
these bursts of dissipation, short-lived and localized, have a measurable impact on molecular abundances in the diffuse medium, and 
how the chemical enrichment they drive compares to observations, are the questions we address here. 

Methods. We further explore a hybrid method to compute the chemical and thermal evolution of a magnetized dissipative structure, 
under the energetic constraints provided by the observed properties of turbulence in the cold neutral medium. For the first time, we 
model a random line of sight by taking into account the relative duration of the bursts with respect to the thermal and chemical 
relaxation timescales of the gas. The key parameter is the turbulent rate of strain a due to the ambient turbulence. With the gas density, 
it controls the size of the dissipative structures, therefore the strength of the burst. It also sets the relative importance of viscous 
dissipation and ion-neutral friction in the gas heating and chemical enrichment. 

Results. For a large range of rates of strain and densities, the models of turbulent dissipation regions (TDR) reproduce the CH + column 
densities observed in the diffuse medium and their correlation with highly excited H2. They do so without producing an excess of CH. 
As a natural consequence, they reproduce the abundance ratios of HCO + /OH and HCO + /H 2 0, and their dynamic range of about one 
order of magnitude observed in diffuse gas. Large C 2 H and CO abundances, also related to those of HCO + , are another outcome of the 
TDR models that compare well with observed values. Neutral carbon exceeds the abundance expected at ionization equilibrium, in 
agreement with fine- structure line observations. The abundances and column densities computed for CN, HCN and HNC are one order 
of magnitude above PDR model predictions, although still significantly smaller than observed values. The dependence of our results 
on the rate of strain and density reveals that the chemical enhancements are in better agreement with observations if the dissipation is 
I dominated by ion-neutral friction, involving shear structures of thickness ~ 100 AU. 

Key words. Astrochemistry - Turbulence - ISM: molecules - ISM: kinematics and dynamics - ISM: structure - ISM: clouds 

o3 ■ 1 ■ Introduction (Miville-Deschenes et al. 2003), 

(3) its spatial and velocity structure is even a greater challenge 

The diffuse medium has a major contribution to the mass of the since il has to reconcile the existence of structures observed at 

interstellar medium (ISM) in galaxies like the Milky Way and a11 scales in emission and a remarkable similarity of line profiles 

as such is a key player in the star formation process. Although observed in absorption (see the discussion in Liszt & Lucas 1998 

it is the first component of the ISM to have been discovered, who question the elusive dimensionality of the diffuse medium), 

and later on extensively analyzed through absorption measure- since man y tracers of the kinematics and small-scale structure of 

ments of atoms, ions and molecules (see the review of Snow & the diffuse ISM are molecular lines, clues are likely to be found 

McCall, 2006), its structure and properties remain a challenge in in its surprisingly rich, but poorly understood, chemistry, 
many respects: 

(1) long thought to consist of two stable phases - the warm and One major puzzle in this chemistry, raised by the detection of 
cold neutral medium (WNM at temperatures T ~ 8000 K and CH + in almost every line of sight sampling the CNM, persisted 
CNM at T ~ 100 K) in thermal pressure equilibrium - a sig- for decades because no formation path was found to be efficient 
nificant fraction of its emission is now detected at temperatures enough in the diffuse ISM (Black, Dalgarno & Oppenheimer 
covering the whole range between those of the CNM and WNM 1975, Black & Dalgarno 1977). It is rooted in the fact that in 
(Heiles & Troland 2003), such diffuse gas, CH + forms via a highly endoenergetic reaction 

(2) the CNM is turbulent with supersonic velocities, yet the C + + H2 (AE/k = 4640 K) unlikely to proceed at the low tem- 
velocity and density power spectra carry the signature of the peratures of the CNM. Similarly, one way to activate the oxygen 
Kolmogorov power spectrum for incompressible turbulence chemistry in the diffuse medium, and therefore the formation of 
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OH and H2O, involves the reaction of O + H2 which has an en- 
ergy barrier of 2980 K. 

A possibly related issue is the existence of rotationally ex- 
cited H2 in the diffuse gas. FUSE observations have revealed 
large populations in the J > 2 levels of H2, inconsistent with 
fluorescence generated by the ambient UV field (Sonnentrucker 
et al. 2003; Gry et al. 2002; Nehme et al. 2008a,b; Martin-Zaidi 
et al. 2005; Lacour et al. 2005). This excited H2 could be local- 
ized in circumstellar materiaQ but it also has been detected in 
the direction of late B stars, devoid of circumstellar matter, as in 
the data of Gry et al. (2002). Stellar UV photons are therefore 
unlikely to contribute significantly to the UV-pumping of H2. In 
particular Lacour et al. (2005) find an increase of the Doppler 
parameter of the H2 lines with J, supporting the existence of a 
warm component that cannot be heated by UV photons. They 
argue that this warm component cannot be due to H2 formation 
pumping, as proposed by Sternberg & Dalgarno (1995) in dense 
PDRs, because it would require an H2 formation rate larger than 
that inferred from observations, and would not reproduce the ob- 
served column densities of CH + found to correlate with excited 
H 2 (Spitzer et al. 1974; Snow 1976, 1977; Frisch & Jura 1980; 
Lambert & Danks 1986). 

ISO-SWS observations further support the possible exis- 
tence of a small fraction of warm gas in the Galactic diffuse 
medium by revealing its pure rotational line emission (Falgarone 
et al. 2005). Interestingly, the ratio NQl^varmlN^ ~ 2 x 10~ 4 , 
where N(H2) wa rm is the H2 column density in levels J > 3, is the 
same across the Galactic diffuse medium as in the direction of 
nearby late B stars. Recent Spitzer observations have confirmed 
the ISO-SWS line flux values (Verstraete et al. in preparation). 

Both the observed abundances of CH + and column densi- 
ties of rotationally excited H2 suggest that large amounts of 
suprathermal energy are deposited in the cold diffuse medium. 
One obvious reservoir of suprathermal energy in the ISM is its 
turbulent kinetic energy. Attempts at incorporating this energy 
in the chemical networks of magneto-hydrodynamical (MHD) 
shocks have been partly successful at reproducing the observed 
properties of the diffuse medium (Pineau des Forets et al. 1986, 
Draine & Katz 1986, Flower & Pineau des Forets 1998). Other 
routes have been explored, involving dynamic interactions of 
the gas and the star cluster in the Pleiades (White 1984, 2003; 
Ritchey et al. 2006), turbulent transport between the WNM 
and CNM (Lesaffre et al. 2007) and turbulent dissipation tak- 
ing place in regions of large velocity shears. Turbulence be- 
ing intermittent in space and time (see the review of Anselmet, 
Antonia & Danaila 2001), velocity shears may locally be large 
enough to drive large local heating rates and trigger the endoen- 
ergetic reactions of carbon and oxygen chemistries in the diffuse 
ISM (Falgarone, Pineau des Forets & Roueff 1995). Along these 
lines, Joulain et al. (1998, hereafter J98) have explored the role 
of ion-neutral decoupling induced, in the weakly ionized diffuse 
medium, by the sharp gas accelerations in the regions of largest 
velocity-shear and its impact on ion-neutral chemistry, in partic- 
ular the formation of CH + and HCO + . Falgarone et al. (2006) 
have analysed the thermal and chemical relaxation phase in the 
evolution of a gas chemically enriched in a dissipation burst. 

The observational data probing the molecular diversity and 
richness of the diffuse ISM are not restricted to CH + and HCO + . 



1 Here, we deliberately overlook most of the early Copernicus results 
on H2 high-./ lines obtained in the direction of hot stars (e.g. Spitzer 
et al. 1973; Savage et al. 1977; Shull & Beckwith 1982) that led to 
the conclusion that excited H2 absorption is occurring in circumstellar 
material heated/shocked by the star itself. 



In their long-lasting effort dedicated to unravelling molecular 
abundances in the diffuse ISM, Liszt & Lucas (2002 and refer- 
ences therein) have provided us with invaluable constraints. Not 
only did they show that the abundances of several molecules stay 
proportional to each other, with very well defined abundance ra- 
tios, but they found that the column densities of these molecules 
vary by more than one order of magnitude across clouds that all 
have about the same total hydrogen column density, correspond- 
ing to diffuse and translucent clouds. They also revealed the re- 
markable similarity of HCO + and OH line profiles, all the more 
surprising for an ion and a neutral species, differently coupled 
to the magnetic field. On the contrary, the high-spectral resolu- 
tion spectra of Crane, Lambert and Sheffer (1995) convincingly 
showed that the CH + line profiles are definitely broader and less 
Gaussian than those of CH, along the same lines of sight, while 
Lambert, Sheffer & Crane (1990) found that, in the direction of 
£ Oph, the CH line profiles could be seen as the sum of a broad 
component similar to the CH + profile, and a narrow one, close 
to that of the CN line. A similar result was obtained by Pan et al. 
(2004, 2005) towards stars of the CepOB2 and CepOB3 regions. 
These sets of results suggest that the velocity field is involved 
in the origin and the evolution of these molecules, and does so 
differently for each species. 

It is therefore challenging to compare these available obser- 
vations with models of a random line of sight across the dif- 
fuse ISM, where active dissipation bursts coexist with others in 
their relaxation phase. In particular, the possibility that a number 
of transient events may dominate the observed molecular col- 
umn densities has never been addressed. This is what we do in 
the present paper. We restrict our study to densities lower than 
«h = 200 crrT 3 because, as will be seen, at higher densities and 
for the turbulent energy observed in the CNM, turbulent dissi- 
pation does not heat the gas enough to open the endoenergetic 
barriers mentioned above. We cannot rule out, though, rare dis- 
sipation bursts of exceptional intensity that would be able to heat 
still denser gas to the required temperature. We extend the previ- 
ous studies in a way that allows us to explore the parameter do- 
main, in particular those characterizing the ambient turbulence. 
We also extend the chemical network. We model a random line 
of sight across the diffuse medium and compare the predicted 
column densities of a variety of molecular species to the obser- 
vations. The dynamic steady state is computed and described in 
Section 2, the chemistry in the active dissipation and relaxation 
phases is presented in Sections 3 and 4. The modelling of a line 
of sight is discussed in Section 5. Comparisons of computed col- 
umn densities with observations, as well as excitation diagrams 
of H2, are shown in Section 6 and these results are discussed in 
Sect. 13 

2. Steady state of a magnetized vortex in a 
weakly-ionized diffuse gas 

2.1. The neutral flow 

Turbulence in the diffuse ISM is supersonic with respect to its 
cold phase, the CNM. Supersonic turbulence dissipates in shocks 
and regions of large velocity shear (Kritsuk et al. 2007). Their 
respective importance has been studied in numerical simula- 
tions. Porter et al. (1992; 1994) and Pavlovski, Smith & Mac 
Low (2006) showed that most of the turbulent kinetic energy 
is rapidly transferred to high wavenumber non-compressible 
modes, once the shocks generated by supersonic turbulence have 
started to interact, reducing the role of compressible modes 
(shocks) in turbulent dissipation. In the so-called quiescent ISM 
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(i.e. far from star forming regions), the smooth observed line- 
shapes support the view of a turbulence devoid of strong shocks 
(Falgarone et al. 1994) implying that the dissipation preferen- 
tially occurs in shear-layers. 

Dissipative structures are modelled as shear-layers belonging 
to a solution of the Helmholtz equation for vorticity, close to 
the Burgers vortex adopted in J98: this analytical solution has 
the merit that it has only two free parameters that describe the 
balance between the stretching action of the large scales and the 
diffusion of vorticity across the vortex edge, at small scale. It 
provides an analytical framework in which we can compute the 
effect of partial decoupling between ions and neutrals upon the 
steady state configuration of velocity and magnetic field and thus 
explore the effect of both the ion-neutral friction and viscous 
heating upon the chemical network. 

The modified Burgers vortex is an axisymmetric solution 
elaborated in atmospheric sciences by Nolan & Farrell (1998). 
It is identical to the Burgers vortex at small radii r, in cylindrical 
coordinates (r, 6, z), and differs from it at large radii in the sense 
that the radial inflow velocity is not divergent at infinity: 



u r (r) = ~-r ■ 



-fir 1 



(1) 



where a is the turbulent rate of strain (in s _1 ) and /3 (in cm -2 ) 
describes the cut-off of the radial velocity. The axial velocity u z , 
the vorticity a> z and the orthoradial velocity ug are inferred from 
the Helmholtz and continuity equations: 



u z (r) =az-e-^ -(l -/3r 2 ), 



_ _SL l- e -*" 

(jD z (r) — (l>() ■ e p t 



Ug(r) 



i r 

r'u z (r')dr' 

r Jo 



(2) 

(3) 
(4) 



where uq is the peak of vorticity and v is the kinematic viscosity. 
Any vortex is therefore entirely defined by three parameters, a, 
P and (Oq, 

Note that, according to the radial dependence of the vorticity, 
the same equilibrium vortex radius ro as for the Burgers vortex 
can be defined, 



r\ = 4v/a 



(5) 



involving the two quantities that act on it, the rate of strain a and 
the viscosity v. Accordingly, the vortex crossing time 



J r u r (r') a 



Ik) 



(6) 



for any constant k < 1 depends only on the rate of strain, while 
the vortex period, defined as 



ue(ro) 



depends on v, a, and coq. 

Because vorticity is radially non uniform there is a differ- 
ential rotation of the fluid within the structure which induces a 
viscous dissipation rate: 



r =yvi^ 

Lm A* ''ax. 



(8) 



where nij is the stress tensor. This rate is written in cylindrical 
coordinates: 



du e 



uo 
r 



du, 
~dr~ 



+ 2 



du r 
~dr~ 



+ 2 



+ 2 



dz 



(9) 



where rj — pv is the dynamic viscosity in a gas of density p. It 
is computed for hydrogen atoms using the Kay & Laby (1966) 
tables of physical and chemical constants : r\ — 6 x 10~ 6 T^ 2 g 
cirT 1 s _1 where 7j is the gas kinetic temperature. 

In all the following, it is assumed that the fluid description 
of the gas motions is justified because the mean free path of H 
atoms Ah-h = 0.23(«h/50 cirT 3 ) -1 AU in the diffuse medium, 
for a H-H elastic collision cross sectiorfl o"h-h = 5.7 x 10~ 15 
cm 2 (Spitzer 1978), is smaller than all the lengthscales involved 
in the model. 



2.2. Interstellar constraints on the vortex parameters 

As said above, each vortex is defined by a set of three inde- 
pendent parameters, a, /3 and wq. The cut-off parameter j3 has 
no influence either on the dynamics or on the chemistry of the 
structure as long as u r is small compared to ug >max . Hereafter, f3 
is chosen in order to satisfy this condition. In other words, tur- 
bulent dissipation in the vortex happens through radial vorticity 
distribution, not through a radial flux of matter. In these condi- 
tions, the dominant contribution to the viscous dissipation is the 
first term in Eq. (0. 

Numerical simulations of incompressible turbulence 
(Jimenez 1997) and experiments (Belin et al. 1996) have shown 
that the maximum tangential velocity of filaments of vorticity 
UQ,max ~ wq7"o * s °f me order of the rms velocity dispersion of 
the ambient turbulence cr,„^. Since the equilibrium radius ro 
is set by the rate of strain a (Eq. |5), is also determined. In 
the case of interstellar turbulence, it implies that the orthoradial 
velocity in the vortex is supersonic with respect to the cold gas. 
It is noteworthy that slightly supersonic Burgers vortices have 
been found in experiments of rotating magnetized plasmas by 
Nagaoka et al. (2002) in conjunction with an inner density hole. 
Moreover, as will be seen later, the gas being violently and 
rapidly heated in the layers of largest orthoradial velocity, the 
Mach number there is likely to drop below unity. The only free 
parameter left in the vortex description is therefore the rate of 
strain a, although the gas density is in fact a free parameter that 
determines the vortex size, through the kinematic viscosity v 
(Eq.0. 

2.3. Magnetic field configuration and ionized flow 

The configuration of the magnetic field and ions reached once 
the above vortex has developed in the partially ionized gas is 
numerically computed. The ions are initially at rest, threaded by 
a uniform magnetic field parallel to the z axis, Bo = fiok. The 
ions are predominantly C + , the neutrals are mainly composed of 
(7) H and H2 and the ionization degree, x — 2 x 10~ 4 , is weak0. At 
t — 0, the ions are suddenly put into motion by friction with the 
vortex that developed in the neutral gas i.e. the three components 
of the neutral gas velocity u„ are those of the vortex given by 



For comparison the H2-H9 elastic collision cross section is 
o-h 2 -h 2 ~ 3 1(T 15 cm 2 (Monchick et al. 1980). 

3 Computed for a diffuse molecular cloud of density «h ~ 50 cm 3 , 
temperature 7* ~ 100 K, illuminated by the standard interstellar radia- 
tion field (ISRF). 
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Eqs. (1), (2) and (4). Boundary conditions are provided by the 
assumption that the vortex has a finite length Ly, apodised over 
a length Cy- 

The alignment of to with the ambient magnetic field is sup- 
ported by the results of numerical simulations. Brandenburg et 
al. (1996) showed that in MHD turbulence, magnetic field and 
vorticity vectors tend to align with each other. More recently 
Mininni et al. (2006a,b) observed a similar behaviour in their 
1536 3 numerical simulations. 

Under these assumptions, we compute the two-dimensional 
time-dependent evolution of the ion velocity Ui and the mag- 
netic field B. We neglect the retro-action of the ions upon the 
neutral motions because, for densities in the range 10-200 cirT 3 
and an ion-neutral drift velocity comparable to u„ (see Fig. 
1), the friction force Fj n they exert on the neutrals is negli- 
gible compared to the advection force in the vortical motion: 
Fj„ ~ 10~ 3 (//10AU) x p„ u n .Vu n , I being the spatial scale for the 
variation of u„, in the range 10 to 100 AU. The neutrals velocity 
components are therefore those of the vortex (Eqs. (1), (2) and 
(4)) at any time. 

In the interstellar medium B is frozen in the charged fluid 
(Spitzer 1978) and its evolution is simply written: 



SB 

— + V x (B x u0 = 0. 

ot 



(10) 



Neglecting the pressure gradients in the evolution equation 
of the ionized flow (this assumption is justified in Sect. |7]i leads 
to: 



<9ui (cry), 
— +(ui-V) Ui = - — — 

at (jj„ + no 



p n (u„ -Ui) + -^(VXB)X B(l 1) 



where p n and ju, are the mean mass per particle of the neutrals (H, 
H2) and ions (mostly C + ) respectively. {crv) in = 2.2 x I0~ 9 cm 3 
s _1 is the momentum transfer rate coefficient between the ion- 
ized and neutral fluids calculated by Flower & Pineau des Forets 
(1995, Appendix A), and close to the Langevin rate. 









M 


Mi 


M 2 


M 3 


Magnetic field 


B 




10 


5 


10 


10 


Vortex length 




AU 


200 


200 


100 


200 


Apodisation length 


C v 


AU 


100 


100 


50 


50 



Table 1. Parameters of the four models of Fig. [T] The density 
« H and the turbulent rate of strain are fixed: n H = 50 crrT 3 , a = 
V 



We integrate Eq. dTOb and Eq. (fTTT > by means of a two dimen- 
sional implicit scheme using the Alternating Direction Implicit 
method (ADI). To validate our approach we also use two other 
integration schemes: an explicit and an implicit without the ADI 
method. The results of our 300 x 200 points grid simulations are 
displayed in Fig.Q] Panel (a) shows the evolution of |m,j_|„ mv , the 
maximum ion velocity component perpendicular to the axis z, 
as a function of time. The curves correspond to the four models 
presented in TableQ] Panels (b) and (c) show the orthoradial and 
axial components of B for the model Mo at t — 100 yr. 

Fig- Eh shows that the ions, initially at rest, are put into mo- 
tion by the friction force from the neutrals in the vortex. This 
motion (including its boundary conditions) generates an ortho- 
radial component Bg of B and a gradient of B z . These terms in 
turn induce magnetic tension and pressure gradient (see Eq.fTOb. 
two forces resisting the orthoradial entrainment exerted by the 




10 ™ 

20 30 
r (AU) 40 



z (AU) 



Figure 1. Magnetic properties of the vortex. Panel (a): evolution 
of \iij ± \ max , the maximum ion velocity component perpendicular 
to the axis z, as a function of time. The different curves corre- 
spond to the models M , Mi, M 2 and M 3 , (see Table [TJ. Panels 
(b) and (c): orthoradial and axial components of the magnetic 
field at t = 100 yr for the model Mo as functions of r and z- 



neutrals. After ~ 100 yr, |Ujj_| mav < 0.4 km s which is small 
compared to \u n± \ in the vortex (see Fig. [2] in Sect. 13.21 1 for all 
models. A steady state is reached in which the ions are almost 
back to rest and the magnetic field slightly helical (Fig. QJ)). A 
large and steady state drift is set between the ion and neutral 
orthoradial velocities with an amplitude close to the orthoradial 
velocity of the neutrals in the vortex. 

Such a drift has a deep impact on the chemistry of the gas, as 
was shown by J98, and contributes to the dissipation of its tur- 
bulent energy. The additional heating term due to the ion-neutral 
friction is written: 



PnPi 



■ (oy)i, 



Mi 



■ (Ui - u„) 



(12) 



p„ + [li jJln + Mi 

We also found that there is only a very small ion-neutral drift 
in the z direction because the magnetic field is only slightly he- 
lical. Since in addition the contributions of the spatial variations 
of u z to the heating term T nn are negligible (Eq. |5J, we consider 
the vortex as invariant along the axis z and restrict our study of 
the spatial and time dependences to those occurring radially. 

In the next sections we focus on the rapid thermal and chem- 
ical evolution of the gas trapped in such a steady state structure, 
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and follow its thermal and chemical relaxation, once the vortex 
has blown-up. 



—r- = -r ^riikTi \ = Bi + B e + —T in + — r e 
at at \1 J fij fi e 



(14) 



3. The active phase 

3. 1 . Numerical modelling 

As in J98, we follow the Lagrangian evolution of a fluid parti- 
cle trapped in the steady state vortex configuration. Because the 
vortex crossing time r c (see Sect. [XT) is comparable to the chem- 
ical timescales, we compute non-equilibrium chemistry coupled 
to the time-dependent thermal evolution. The initial gas temper- 
atures and molecular abundances are those of a steady state dif- 
fuse cloud of density «h, illuminated by the ambient interstellar 
radiation field (ISRF) (Draine 1978) scaled by the factory, and 
shielded by the extinction A v . The cosmic ray ionization rate £ 
and the elemental abundances are given in Table [2] 

The neutrals and the ions are treated separateljQ and a fluid 
particle is defined at each time by its position r, the neutral and 
ionized velocity fields u n and Ui, the mass densities p„ and p,-, 
the temperatures of the neutrals T„ and ions 7f, and the abun- 
dances n(X) of each species. The system therefore comprises 1 1 
dynamic time-dependent variables (r, n,-, n„, pi, p„, r,-, T„, u nr , 
u„g, Ui r , u#)). Our chemical network originates from the Meudon 
PDR code (Le Petit et al. 2006). It incorporates 90 species inter- 
acting through 1524 reactions. Those include the formation of 
H2 on dust, the photoprocesses and the processes induced by the 
cosmic rays. We also compute the time-dependent evolution of 
the populations of the 18 first ro-vibrational levels of H2 (corre- 
sponding to T ex < 10 4 K). 



Density 


n H crrT J 10-200 




Radiation field 


X 


1 




Extinction 


A v mag 0.1-1 




Cosmic ray ionization rate 


( s- 


1 3xlO~ 17 


a 


Elemental abundances 




[X] / [H] 




Helium 


[He] 


1.00 x 10" 


-] 


Carbon 


[C] 


1.38 x 10" 


-4 


Oxygen 


[O] 


3.02 x 10" 


-4 


Nitrogen 


[NJ 


7.94 x 10" 


■5 


Sulfur 


[S] 


1.86 x 10- 


-5 


Iron 


[Fe] 


1.50 x 10- 


-8 



Table 2. Physical conditions and elemental abundances of the 
gas in the TDR models, a - Dalgarno (2006). 



The system of variables is therefore driven by a set of 
119 first-order coupled differential equations that are integrated 
along the fluid particle trajectory. To ensure that the time step is 
consistent with the variations of all dependent variables we use 
the DVODE differential equation solver (Brown et al. 1989). 

The evolution of the thermal energy densities U„ and t/, is 
given by: 



dU„ 
dt 



= ~ ||nnWJ = b„ + r,„, + 



en 



(13) 



4 The ions and the electrons are treated as a unique fluid at a temper- 
ature Tj because: (1) the ion-electron velocity drift induced by the mag- 
netic field fluctuations in the model is ~ 1 cm s -1 , negligible compared 
to the thermal velocities, (2) the ion-electron temperature equipartition 
time is ~ 0.1 (r/1000K) 3/2 (n,/3 lO^cm" 3 )- 1 yr (Spitzer 1978), also 
negligible compared to the dynamic timescales. 



where B„, Bi and B e are the sums of all the heating and cool- 
ing rates of the neutrals, the ions and the electrons respectively, 
not induced by turbulent dissipation, and Y„„ and F,,, the heat- 
ing rates induced by turbulent dissipation defined in the previous 
section. Y en is the heating rate due to the electron-neutral drift, a 
term taken into account in the code but negligible compared to 
T nn and r,„. The cooling rates include the radiative de-excitation 
of the ro-vibrational levels of H2, of the fine structure levels of 
C + , C and O and of the rotational levels of H 2 0, OH and CO. 

Lastly, since ions and neutrals are decoupled, we use the ap- 
proximation of Flower et al. (1985) for the calculation of the 
chemical rate coefficients. The cross-section of a 2-species reac- 
tion is integrated over a Maxwellian velocity distribution at an 
effective temperature 



miT 2 + m 2 T i 1 



'«1'«2 2 



mi + nil 



3k mi + nil 



(15) 



where ni\, ni2, T\ and Tz are respectively the masses and the tem- 
peratures of the 2 reactants and u D their relative drift velocity. 



3.2. Thermal evolution of the gas 

Fig. [2] displays the main properties of a reference model where 
a — 3 x 10~ n s , n H = 30 cirT 3 and Ay = 0.4 mag. The vortex 
has an equilibrium radius ro = 38 AU and generates an average 



turbulent heating rate, T tur b 
as: 



= 3.4 x 10~ 23 er 



g cm 



turb 



2/(Kr Q ) 



f 

Jo 



[T nn (r) + T in (r)]rdr. 



s 1 defined 



(16) 



The integration domain extends to the radius Kro where the tur- 
bulent heating has no significant influence on the gas tempera- 
ture and chemistry {K ~ 5). 

In the model presented here, the heating rate is dominated 
everywhere by ion-neutral friction (Fig. |2^). For higher values 
of the turbulent rate of strain a, ro decreases (see Sect. 15.3b and 
the importance of the viscous dissipation increases because u„g 
is fixed by the ambient turbulence. 

The gas in the vortex never reaches thermal balance and the 
thermal inertia is visible by comparing Figs.|2^ and|2}; where the 
peak temperature of the fluid particle is reached a few hundred 
years after the peak of the heating rate. Emission in the pure 
rotational lines of H2 is by far the dominant coolant in the layers 
where T„ > 200K (Fig. [2};) while the cooling rate due to the 
ionized carbon C + decreases in the warmest layers. This is due 
to the chemical evolution (see Sect. |3.3I >. 

Last, some chemical clues are provided in Fig. [2};. The 
neutral-neutral reactions only depend on the temperature T n 
while the ion-neutral reactions depend on the ion-neutral drift. In 
Eq. Sl5[ , the second term of the right hand side reaches 1000K 
as ud ~ 3 km s~', an effective temperature higher than the peak 
kinetic temperature in the vortex. A comparison of the shapes of 
the orthoradial velocity u„g and the neutral temperature T„ shows 
that the endo-energetic ion-neutral chemistry is activated earlier 
in the fluid particle evolution than the neutral-neutral chemistry. 
For each type of vortex, the relative importance of those two 
chemical regimes is different. 
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Figure 2. Vortex physical properties as functions of the radius 
(bottom axis) or time (top axis, arbitrary origin) for the refer- 
ence model: a — 3 x 10~ n s _1 , n H = 30 cirT 3 and A v - 0.4 mag. 
Panel (a): The heating terms T nn due to the viscous dissipation, 
Tin due to the ion-neutral friction and T pn - e due to the photo- 
electric effect. Panel (b): The dominant cooling terms due to the 
radiative desexcitation of H2 (pure rotational lines), C, O and C + 
(fine-structure lines). Panel (c): The temperature and orthoradial 
velocity of the neutrals. 



3.3. Chemical evolution of the gas 

The chemical evolution of the gas during the vortex active phase 
is similar to that reported in J98, although the chemical network 
is updated and includes nitrogen- and sulfur-bearing molecules. 
The outline of this network is given in AppendixIClwhere we dis- 
play the main production and destruction routes of the molecules 
of interest (1) in the ambient diffuse medium («h = 30 cm 4 and 
Ay - 0.4 mag) and (2) in the vortex for the reference model at a 
radius r — ro. 

The most important reaction route opened by the dissipative 
structure is the endothermic hydrogenation of C + : 

C + +H 2 ^CH + + H - AE/k = 4640/:. (17) 

Besides the direct production of CH + , this reaction is respon- 
sible for most of the chemical richness of the vortex as shown in 
Fig. |C. 21 (Appendix [Ct: it enhances the production of CHt via 
the successive hydrogenation by H2 of CH + and CHJ. CH3, in 
turn, enhances the production of molecules including CH and 



C via the dissociative recombination with electrons, HCO + and 
CO via 

CH+ + O -> HCO + + H 2 (18) 

and CN, HCN and HNC via 

CH+ + N -> HCN + + H 2 . (19) 

The production of CH3 is also at the origin of C 2 H and CS since 
these molecules are both products of CH (through the reactions 
CH + C + — > Cj + H and CH + S + -» CS + + H respectively). 

The second main reaction, absent in the ambient medium, 
which plays an important role in the chemical evolution of the 
vortex is: 

O + H 2 -> OH + H - AE/k = 2980K. (20) 

Besides the direct production of OH it triggers the production of 
H 2 via 

OH + H 2 ->H 2 + H -AE/k= 1490K (21) 
and 2 via 

OH + O -> 2 + H. (22) 

FigS displays the evolution of several relative abundances 
in the magnetized vortex. The impact of the turbulent heating 
is such that most species abundances rise from 2 to 5 orders of 
magnitude within the structure. The formation of vortices in the 
turbulent gas flow therefore has specific chemical signatures that 
we expect to observe in the diffuse medium. One remarkable 
example is HCO + . Since this molecule is a direct product of CHt 
(in the vortex) it becomes related to almost all the species whose 
abundance is enhanced in the vortex. 

Time (and space) stratification is also visible in Fig[3j 
the chemical enrichments occur successively, because the ion- 
neutral chemistry is triggered earlier than the neutral-neutral 
chemistry (see previous section) and because the chemical in- 
ertia of each species is different. 

4. The relaxation phase 

4.1. Numerical modelling 

Once the burst of turbulent dissipation is over, some chemical 
signatures imprinted in the gas persist for several thousand years 
as shown by Falgarone et al. (2006). To compute the chemical 
and thermal relaxation of the gas, we assume that, once the vor- 
tex has vanished, the gas is dynamically frozen: u = 0. The 
previous Lagrangian approach is switched to Eulerian, and we 
compute the time-dependent evolution of each cell in the vortex. 
The initial conditions of the relaxation are the conditions of the 
active stage at every position. The thermal equations ( fT3] > and 
(fl4l are still valid with T„„ = and F,„ = 0. 

While the numerical code has been conceived to treat iso- 
baric or isochoric relaxation, all the results presented in this pa- 
per were obtained assuming isochoric relaxation, because it al- 
lows us to better disentangle what is due to the chemistry itself 
from what is due to the gas density. In particular, it shows more 
clearly the role of the relaxation timescales of the molecules, i.e. 
only driven by the chemical network and the thermal evolution, 
independently of the gas density. 
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Figure 3. Fractional abundances relative to «h = «(H) + 2w(H2) 
of selected species as functions of the radius (bottom axis) or 
time (top axis, arbitrary origin) for the reference model: a = 
3 x KT 11 s -1 , h h = 30ctrT 3 and Ay = 0.4 mag. 



Figure 4. Vortex thermal features during the relaxation phase, 
as functions of the radius and time (after the vortex blow-up) 
for the reference model: a = 3 x 10~ u s _1 , «h = 30 ctrT 3 and 
Ay = 0.4 mag. An isochoric relaxation is assumed. Top panel: 
the main cooling terms, i.e. the radiative desexcitation of H2 and 
C + . Bottom panel: the temperature of the neutrals. 



4.2. Thermal evolution of the gas 

Fig. E] displays the main cooling rates (top panel) and gas tem- 
perature (bottom panel) as functions of time and position in the 
vortex (after the vortex blow-up). It shows that as in the active 
phase, the cooling rate during the relaxation phase is dominated 
by the emission in the pure rotational lines of H2. In the model 
presented here («h = 30 cm 4 ), other cooling agents (mainly 
C + ) become dominant at t ~ 10 4 yr. 

4.3. Chemical evolution of the gas 

During relaxation, the cooling of the gas causes all the endo- 
energetic reactions triggered during the dissipative burst to slow 
down and stop, one after the other. The gas loses its chemical en- 
richment at a speed that depends on the molecular species. This 
is illustrated in Fig.|5]that displays the time-dependent evolution 
of the column density Nvr(X, t) of selected species X integrated 
across the vortex. 

For most species (CO, C2, OH, H2O, C2H...) the signature 
of the turbulent dissipation persists over more than 10 3 yr. The 
characteristic timescales (e-folding times) differ between species 
by more than a factor of 30 (~ 2 x 10 2 yr for CH + and 7 x 10 3 yr 



for CO). The existence of the relaxation phases modifies the cor- 
relations between molecular abundances. One example is pro- 
vided by HCN and C2H that follow a similar enhancement in the 
active phase (see Fig. [3} and have markedly different behaviours 
in the relaxation phase. 

The richest phases are not necessarily those contributing 
the most to the observable column densities because they have 
short lifetimes. In the next section, we detail how we take time- 
variability into account in our modelling of a random line of 
sight across the diffuse medium. 

5. Modelling of a line of sight 

In turbulent flows, the spatial distribution of the regions of high- 
est dissipation rate (extrema of velocity shear, extrema of nega- 
tive velocity divergence i.e. shocks) is far from space-filling, one 
of the aspects of the intermittent nature of turbulence. The filling 
factor of these extrema has been computed in numerical simula- 
tions of mildly compressible (Pety & Falgarone 2000) or super- 
sonic MHD turbulence (Pavlovski et al. 2006, Pan & Padoan 
2008). More recently, Moisy & Jimenez (2004) have shown that 
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Figure 5. Column densities of selected species, integrated across 
the reference vortex, as functions of time during an isochoric 
relaxation phase. 



the regions of intense vorticity tend to form filaments, while re- 
gions of most intense dissipation rather form sheets or ribbons, 
all of them being organized in clusters, probing the organization 
of small-scale intermittent structures. 

For our purpose, we assume that any line of sight intercepts a 
number of vortices, either active or in their relaxation phase. The 
total number of vortices per line of sight is constrained by the av- 
erage transfer rate of turbulent energy per unit volume available 
in the cascade. 

5.1. Method 

Any line of sight samples three kinds of diffuse gas: (1) the ambi- 
ent medium in which the chemistry is computed as steady state 
UV-dominated chemistry, (2) a number Nva of active vortices 
and (3) a number of relaxation phases related to Nva- We also 
consider that the line of sight is homogeneous and characterized 
by its density n H and uniform shielding A v from the ISRF. Such 
an approximation is useful to test the importance of each param- 
eter (a, «h and Ay) on the final chemical state of the gas. Last, 
we assume that all the vortices in the line of sight are identical: 
they all have the same turbulent rate of strain a. In the following 
section we show why a more realistic description, with a distri- 
bution of rate of strain values, would not provide very different 



results, namely because, as it will be seen, the results depend 
weakly on a. 

The contribution of one active phase of duration Ty to the 
total column density of a species X is Nva(X). The contribution 
of the relaxation phase of that species is computed by assuming 
that the longer its relaxation timescale, the greater the contribu- 
tion of the relaxation phase in the observed column density so 
that the chemical composition of a line of sight is entirely deter- 
mined by the number of active vortices Nva and their lifetime 
Ty. The resulting column density of a species X is: 



N(X) = N VA (X) 



f 

Jo 



N va (X) + 1/t v N VR (X,t)dt 



N M (X) (23) 



where Nm(X) is the contribution of the ambient medium. 

5.2. Constraints provided by the energy available in the 
turbulent cascade 

5.2.1. The number of active vortices 

The number of active vortices Nva in a line of sight and their 
lifetime t v , are constrained by the turbulent energy available in 
the cascade and its transfer rate. In turbulent flows, the transfer 
rate of kinetic energy at scale I of characteristic velocity u\ is, on 
average: 



1 u, 
Sl = 2 P T 



(24) 



The Kolmogorov scaling of turbulence, valid for incompress- 
ible turbulence, postulates that this transfer rate is scale- 
independent. In the highly compressible interstellar medium, 
one would expect this quantity to differ from one scale to an- 
other. Measurements of the internal velocity dispersion of clouds 
of size I combined with their density provide an estimate of 
the turbulent energy transfer rate at this scale. A compilation of 
CO(l-O) line observations of interstellar structures of size rang- 
ing between 10 -2 and 10 3 pc shows that e/ is remarkably in- 
dependent of the scale in the Galaxy and that there is a large 
scatter (by a factor of 100) about a well defined average value 
s obs ~2x 10~ 25 ergcrrrV (Falgarone, 1998; Hily-Blant et al. 
2008). A similar value holds for turbulence in the HI CNM and 
in non-star-forming dense cores (Falgarone 1999). The unifor- 
mity of this value across the local ISM suggests that the turbulent 
cascade encompasses the different regions of the cold medium, 
and that the transfer is driven by turbulence at the same rate, in 
all media, whatever the gas density. 

We adopt the above value of s b s in our model as representa- 
tive of the turbulent energy transfer rate through scales. We thus 
impose that, at any time, the dissipation rate in all the active vor- 
tices in a line of sight is equal to the average energy transfer rate 
in the turbulent cascade, so that: 



Sobs = NvA^turb2KrolL 



(25) 



where L is the depth of the line of sight, inferred from L = 
Nu/nu- This fixes the number of active vortices in a given line 
of sighQ This number and many of the results are therefore pro- 
portional tO £ bs- 



5 Nva also depends on the angle of inclination of the vortices along 
the line of sight. Because the dynamics and the chemistry in a vortex 
do not depend on the axial coordinates z and because the model is ax- 
isymmetric, this angle has no influence on the final results. Nva is thus 
defined for vortices perpendicular to the line of sight. 
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Figure 6. Main physical properties of the TDR models as functions of the turbulent rate of strain a for several values of the density 
«h- Top panels: the vortex heating terms F„„ due to the viscous dissipation (left) and r,„ due to the ion-neutral friction (right). 
Bottom panel: the vortex lifetime Ty. 



5.2.2. The velocity dispersion of the ambient turbulence 



The energy transfer rate depends on the density, velocity dis- 
persion and timescale. We thus need to know the amount of 
turbulent energy available in the CNM component of the dif- 
fuse medium, or the rms turbulent velocity in the CNM. It is 
this quantity that sets the angular velocity of the vortex (see 
Section 2.2). This quantity is difficult to determine on observa- 
tional grounds, because of the mixture of WNM and CNM in 
the HI emission spectra and because of the lack of spatial in- 
formation for the absorption spectra (dominated by the CNM). 
We adopted a rms velocity dispersion cr, llr i, = 3.5 km s _1 for the 
CNM turbulence derived from the HI maps of a high latitude cir- 
rus in the Polaris Flare (Joncas et al. 1992; Miville-Deschenes 
et al. 2003) in which HI emission is well correlated with the 
100;um emission of dust, probing column densities of gas repre- 
sentative of the CNM. This value is consistent with those quoted 
in Crovisier (1981) for the CNM. It is comparable to the geo- 
metric mean of the two smallest values, FWHM = 4.9 and 12.0 
km s _1 inferred by Haud & Kalberla (2007) from the Gaussian 
decomposition of the HI profiles of the Leiden/ Argentine/Bonn 
survey of galactic HI (Hartmann & Burton 1997). Last, this rms 
velocity dispersion is consistent with the approximate equiparti- 
tion between magnetic and turbulent energy inferred by Heiles 
& Troland (2005) from the median value of the magnetic field 
estimated in the diffuse medium B — 6 //G. 



5.2.3. The lifetime of the active phase 

The lifetime of an active vortex, Ty, is controlled by the large- 
scale motions of the ambient turbulence that feed energy into 
the small-scale structures (Moffatt, Kida & Ohkitani 1994). It 
may exceed the period of the vortex as suggested by a variety of 
experiments in incompressible turbulence (e.g. Douady, Couder, 
Brachet 1991). This lifetime is uncertain, though, and its value 
in our model is constrained by energetic considerations. 

We assume that, for simplicity, all the vortices explored in 
the models dissipate the same total energy over their lifetime Ty 

E=n(Kr ) 2 f turb L V Tv (26) 

This constraint fixes the lifetime Ty of the vortex i.e. the time 
during which turbulent dissipation is active. In order to stay in 
the vortex framework we impose that t v is larger than the vortex 
period P, which sets a lower limit to E. The influence of this 
parameter on the results is discussed in AppendixlBl 

5.3. The role of the rate of strain and density under the 
energetic constraints 

Once s„bs, o~turb an d E are given, a line of sight is therefore de- 
fined by only three independent parameters: the turbulent rate of 
strain a, the gas density «h an d the shielding Ay. 

To help the reader understand the chemical results presented 
in the next section, we discuss the roles of the rate of strain and 
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Figure 7. Main physical properties of the TDR models as functions of the turbulent rate of strain a and of the density «h of the gas. 
Top panels: the vortex turbulent heating rate T tu ,b (left) and the vortex maximum temperature T max (right). Bottom panels: the vortex 
radius rrj (left) and the number of active vortices Nva along a line of sight sampling one magnitude of gas (right). 



gas density because the above constraints on the energy dissipa- 
tion rate actually couple a and «h that should be independent pa- 
rameters. This is so because the two energy constraints (transfer 
rate and energy) involve the equilibrium radius ro that primarily 
depends on a but also on n^, via the density dependence of the 
kinematic viscosity (see Section 2.1) These trends are illustrated 
in Figs. |6]and|2] that also display the dependence of several key 
quantities on a and «h- 

Fig. [6] shows that, as expected, F„„ ~ a«h is almost propor- 
tional to a because the higher the rate of strain, the smaller the 
equilibrium radius and the larger the velocity shear. F,„ is al- 
most independent of a and increases with density as F,„ oc m d «^ 
so that, depending on a and bh, two regimes exist: one at low 
density and high rate of strain where the turbulent heating is 
dominated by viscous dissipation, and the other (small a, high 
density) where this heating is dominated by the ion-neutral fric- 
tion. The rate of strain a therefore plays an important role in the 
nature of the warm chemistry triggered in the vortex. This figure 
also displays the run of -ry with the rate of strain for different 
densities, as a result of the constraint on the total energy E dis- 
sipated in each vortex. 

Fig(7]shows that ro is small in the former regime and reaches 
values of the order of 100 AU in the regime where ion-neutral 
friction dominates. The peak gas temperature reached in the vor- 
tex, T max , is also shown: it is higher in the regime where viscous 



dissipation dominates because the orthoradial velocity is fixed 
and thin vortices induce large velocity shears, thus large viscous 
heating. However, the thermal inertia of the gas prevents it from 
reaching much higher temperatures, because the most efficient 
vortices (large a) are short-lived. The range of a and «h explored 
in our study covers these two regimes and we quantify the chem- 
ical effects of turbulent dissipation as it changes from dominated 
by ion-neutral friction to dominated by viscous dissipation. 

Last, we find that the number of active vortices Nva (com- 
puted for a line of sight sampling one magnitude of gas) is 
roughly independent of a and decreases almost as « H 2 as the den- 
sity increases. This is because s obs is fixed (see previous section), 
and because of the combined dependences of ro, T tur f, and L on 
a and n H ( see Eq.|25l>. We note that Nva reaches large values at 
low density (up to several hundreds along lines of sight of sev- 
eral tens of parsecs). However the filling factor of the vortices 



f, = M 



2Kr n 



VA - 



(27) 



never exceeds f v — 4 x 10 , its lowest values being f v ~ 10 at 
high densities and rates of strain. 
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Table 3. Physical and chemical characteristics of two TDR-models defined by their parameters «h, Ay and a. /m(X), /va(X) and 
/vr(X) are the contributions of the ambient medium, the active and the relaxation phases respectively, to the column density N(X) of 
the species X. All models are computed for Ah = 1.8 x 10 21 cirT 2 . The relative abundances are given in the last columns. Numbers 
in parenthesis are powers of 10. 



Parameters 









30 cm" 3 










100 cnT 


i 






A v 






0.1 mag 










0.1 mag 








a 






3 (-11) s" 


i 








3 (-11) s- 


i 






Physical properties 






9.3 (-24) erg.cm^.s" 1 






6.9 (-23)erg.cnT 3 .s-' 






1 amb 






114 K 










58 K 








T 

* max 






1000 K 










807 K 








TV 

ro 






1070 yr 
40 AU 










506 yr 
21 AU 








Ny A 






215 










17 








fr 






2.2 (-2) 










2.9 (-3) 








Chemical properties 


Species 


/iffl 

% 


.MX) 

% 


/vr(X) 
% 


N(X) 
cur 2 


N(X)/Nn 


Im(X) 

% 


fvA(X) 

% 


fv R (X) 

% 


N(X) 
cirr 2 


N(X)/N H 


H 


100 


0" 


0" 


7.0 (20) 


3.9 (-01) 


100 


Of 


Of 


2.8 (20) 


1.6 


(-02) 


H 2 


100 


0" 


0" 


5.5 (20) 


3.1 (-01) 


100 


0" 


0" 


7.6 (20) 


4.2 


(-01) 


«3 + 


81 


5 


14 


2.9(13) 


1.6 (-08) 


93 


2 


5 


1.2(13) 


6.7 


(-09) 


C 


19 


65 


35 


4.4(14) 


2.4 (-07) 


68 


23 


9 


5.6(14) 


3.1 


(-07) 


CH 


5 


86 


9 


8.1 (12) 


4.5 (-09) 


22 


69 


9 


7.9(12) 


4.4 


(-09) 


CH + 


0.2 


96 


4 


2.2(13) 


1.2 (-08) 


2 


95 


3 


1.3 (12) 


7.2 


(-10) 


c 2 


0.5 


72 


28 


4.0(11) 


2.2 (-10) 


3 


46 


51 


1.4(12) 


7.8 


(-10) 


C 2 H 


0.1 


76 


24 


6.3(11) 


3.5 (-10) 


0.4 


61 


39 


3.5 (12) 


1.9 


(-09) 


OH 


52 


17 


31 


8.7(13) 


4.8 (-08) 


63 


18 


19 


2.3 (13) 


1.3 


(-08) 


H 2 


57 


11 


32 


1.3 (13) 


7.2 (-09) 


68 


14 


18 


3.6 (12) 


2.0 


(-09) 


H,0 + 


35 


23 


42 


1.3 (13) 


7.2 (-09) 


32 


40 


28 


1.1 (12) 


6.1 


(-10) 


CO 


57 


16 


27 


2.6(13) 


1.4 (-08) 


67 


12 


21 


3.0(13) 


1.7 


(-08) 


HCO+ 


14 


69 


17 


8.4(11) 


4.7 (-10) 


11 


77 


12 


3.7 (11) 


2.1 


(-10) 


o 2 


67 


11 


22 


3.8 (10) 


2.1 (-11) 


80 


10 


10 


4.1 (10) 


2.3 


(-11) 


SH + 


0.1 


98 


2 


6.2(11) 


3.4 (-10) 


0.2 


98 


2 


2.1 (11) 


1.2 


(-10) 


CS 


3 


85 


12 


1.8(09) 


1.0 (-12) 


11 


69 


20 


7.0 (09) 


3.9 


(-12) 


HCS + 


0.6 


93 


7 


1.1 (09) 


6.1 (-13) 


1 


89 


10 


1.9 (09) 


1.1 


(-12) 


CN 


1 


79 


20 


2.3(11) 


1.3 (-10) 


8 


63 


29 


1.6(11) 


8.9 


(-11) 


HCN 


0.5 


88 


12 


3.3(10) 


1.8 (-11) 


3 


81 


16 


2.6 (10) 


1.4 


(-11) 


HNC 


2 


89 


9 


5.9 (09) 


3.3 (-12) 


17 


74 


9 


4.0 (09) 


2.2 


(-13) 



a For the reason given in Sect. l5.4l fi/fi(H) = /vr(H) = /vm(H 2 ) = /vrQIi) = 0, because the densities of H and H 2 are not modified in the vortex. 



5.4. The respective contribution of each chemical regime 

Three different chemical regimes enter our line of sight mod- 
elling, associated with the active phase, the relaxation phase and 
the ambient medium. Their respective contributions to the inte- 
grated column density of a species X (see Eq. [231 , 

f M (X) = N M (X)/N(X) (28) 

f VA (X) = N VA (X)N VA (X)/N(X) (29) 

Jr»co 
N VR (X,t)dt/N(X)T V (30) 
o 

depend on this species and are given in Tables. [3] &|4] 

For numerical reasons the contributions /va(X) and /vr(X) 
are computed as excesses of the density n(X) in the active and 
relaxation phases respectively, above its value in the ambient 
medium. Hence Nm(X) = /m(X) x N(X) is the column density 
of the species X along a line of sight without vortices of total 
column density Ah = 1-8 x 10 21 cirr 2 . 

Tables. [3] & @] show that all the regimes have a significant 
contribution, although most of the selected species are predomi- 
nantly formed in the vortices. 



6. Results: Comparison with the observations 

We have computed 150 different models exploring the parameter 
space as follows: 

- 6 densities (n H = 10, 30, 50, 80, 100, 200 cm 4 ) 

- 5 rates of strain (a = 1, 3, 10, 30, 100 xlO" 11 s" 1 ) 

- 5 UV shieldings (Ay = 0.2, 0.4, 0.6, 0.8, 1). 

The turbulent rate of strain extends over two orders of magni- 
tude, the largest value corresponding to the condition ro » Ann 
for all densities, the smallest being the limit of validity of our 
analytical approach. 

In this Section, we present the results of our line of sight 
models (named Turbulent Dissipation Regions models or TDR 
models) on the same displays as large sets of observational data. 
We also compare these data with results obtained from PDR 
models (two-side illuminated slabs of uniform density (Meudon 
PDR code, Le Petit et al. 2006)), computed with the same con- 
ditions as in the TDR models: ^ = 1, £ = 3 X 10~ 17 s~' and with 
the same chemical network. 

The amount of gas along the line of sight in the PDR and 
TDR models is normalized to A^h = 1.8 X 10 21 cm 2 because, as 
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Table 4. Same as Table [3] 
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Parameters 









30 cnT 3 










100 cnT 


i 




A v 






0.4 mag 










0.4 mag 






a 






3 (-11) s" 


i 








3 (-11) s" 


i 




Physical properties 


^turb 




8.5 (-24) erg.cnT 3 .s-' 






6.2 (-23) erg. cm 3 .s 1 




T amb 






79 K 










42 K 






T 

*■ max 






830 K 










710 K 






Tv 

ro 






1250 yr 
38 AU 










610 yr 
20 AU 






N VA 

fr 






246 
2.4 (-2) 










19 
3.2 (-3) 






Chemical properties 




Affl 


.MX) 


/roCX) 


N(X) 


N(X)/N H 


f M (X) 


/va(X) 


MX) 


N(X) 


N(X)/N H 


Species 


% 


% 


% 


cnT 2 




% 


% 


% 


cirT 2 




H 


100 


0° 


0" 


2.0 (20) 


1.1 (-01) 


100 


0° 


0° 


7.9(19) 


4.4 (-02) 


H 2 


100 


0" 


0" 


8.0 (20) 


4.4 (-01) 


100 


0" 


0" 


8.6 (20) 


4.8 (-01) 


H 3 + 


77 


6 


17 


5.9(13) 


3.3 (-08) 


92 


2 


6 


1.4(13) 


7.8 (-09) 


C 


19 


56 


25 


1.3 (15) 


7.2 (-07) 


76 


13 


11 


1.4(15) 


7.8 (-07) 


CH 


2 


87 


11 


3.6(13) 


2.0 (-08) 


20 


66 


14 


1.5 (13) 


8.3 (-09) 


CH + 


0.2 


96 


4 


1.3(13) 


7.2 (-09) 


3 


93 


4 


6.0(11) 


3.3 (-10) 


c 2 


0.3 


57 


43 


4.2(12) 


2.3 (-09) 


3 


30 


67 


4.2 (12) 


2.3 (-09) 


C 2 H 


0.1 


68 


32 


7.4(12) 


4.1 (-09) 


0.4 


52 


47 


1.0(13) 


5.6 (-09) 


OH 


47 


24 


29 


8.4(13) 


4.7 (-08) 


34 


31 


35 


1.5 (13) 


8.3 (-09) 


H 2 


50 


20 


30 


1.3 (13) 


7.2 (-09) 


39 


27 


34 


2.4(12) 


1.3 (-09) 


H 3 + 


20 


46 


34 


9.9(12) 


5.5 (-09) 


11 


59 


30 


7.0(11) 


3.9 (-09) 


CO 


45 


15 


40 


1.2(14) 


6.7 (-08) 


35 


11 


54 


9.1 (13) 


5.1 (-08) 


HCO+ 


6 


82 


12 


2.2(12) 


1.2 (-09) 


3 


84 


13 


4.9(11) 


2.7 (-10) 


o 2 


64 


16 


20 


6.9 (10) 


3.8 (-11) 


52 


19 


29 


4.4 (10) 


2.4 (-11) 


SH + 


0.02 


97 


3 


1.7 (12) 


9.4 (-10) 


0.03 


96 


4 


4.9(11) 


2.7 (-10) 


CS 


1 


82 


17 


1.9(10) 


1.1 (-11) 


9 


61 


30 


3.0(10) 


1.7 (-11) 


HCS + 


0.2 


92 


8 


9.2 (09) 


5.1 (-12) 


1 


86 


13 


6.1 (09) 


3.4 (-12) 


CN 


2 


67 


31 


8.8(11) 


4.9 (-10) 


10 


46 


44 


4.7 (11) 


2.6 (-10) 


HCN 


1 


84 


15 


1.0(11) 


5.6 (-11) 


3 


75 


2 


5.7(10) 


3.2 (-11) 


HNC 


2 


87 


11 


2.7(10) 


1.5 (-11) 


18 


69 


13 


1.2(10) 


6.7 (-12) 



a Same as Table [3] 



will be shown later, it corresponds to most of the data collected 
in the visible, UV and radio ranges. 

6.1. Ultraviolet and visible observations towards nearby stars 

An important clue for the understanding of the diffuse ISM 
chemistry is provided by the combined observations of CH and 
CH + ( Gredel et al. 1993; Crane et al. 1995 ; Gredel 1997; Pan 
et al. 2004, 2005; Sheffer et al. 2008). Although CH + is known 
to be overabundant, the column densities of CH are in agree- 
ment with the predictions of the PDR models. Because the two 
species are tightly related by the chemistry (see Figs. IC. 11 & IC.2I 
in Appendix0, the difficulty of any modelling is then to under- 
stand the physics which leads to an enhancement of CH + without 
changing the amount of CH. Another important clue lies in the 
correlation, mentioned in the Introduction, between CH + and the 
pure rotational J > 3 levels of H2. 

The CH + , CH and excited H2 observational data are dis- 
played in Fig. [8] with the prediction of several PDR models (left 
panels) and TDR models (right panels), computed for diffuse gas 
of density between 10 and 200 cm 4 . As mentioned above, we 
restrict our analysis to diffuse gas illuminated by the ambient 
ISRF, and have therefore removed all the data corresponding to 
lines of sight toward hot stars. 



The TDR models predict column densities of CH + in good 
agreement with the observations for low densities («h < 100 
cm" 3 ). This density limit corresponds to the observed average 
turbulent energy to which we scale our computations (Eq.EBTi: if 
e bs was larger, this limit would be larger too. Most gratifying is 
the fact that this is achieved without producing an excess of CH. 
It is so because, in the vortex, CH is a product of CH + (its pro- 
duction is therefore enhanced by the turbulent dissipation, see 
Appendix |Cj but its destruction is also enhanced since it pro- 
ceeds through endo-energetic reactions (CH + H, -AE /k=2200 
K and CH + H 2 , -A£//fc=1760 K). Note that in panels (c,d), be- 
cause of the assumed homogeneity of the gas (uniform Ay and 
«h along the line of sight), the column densities computed in the 
TDR models are proportional to A^h, hence the straight lines. 

The dissipation of turbulent energy is also a plausible expla- 
nation for the excitation of the pure rotational levels of molec- 
ular hydrogen: the TDR models reproduce not only the CH + to 
HI(./ = 4) column density ratio but also a significant fraction of 
the large dynamic range (2 orders of magnitude) over which the 
correlation is observed (Fig.[8]e,f). We recall that all the column 
densities of species dominated by the warm chemistry are pro- 
portional to s„b s and may be under- or over-estimated by about a 
factor of 10. 

An encouraging result is also shown in Fig. [8](g,h). Gry et 
al. (2002) and Lacour et al. (2005) found that the H2 excitation 
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diagrams obtained in the direction of stars later than B2 cannot 
be reproduced with the stellar UV excitation only. For compar- 
ison we display a PDR and a TDR model, both computed for 
«h = 100 cm 4 , a value close to the optimal density found by 
Nehme et al. (2008b) toward HD 1 02065 . 

The dependence of the TDR models on the density signif- 
icantly differs from those of the PDR models. They share the 
same dependence of the chemistry on gas density but the TDR 
models have the additional dependence on density of the number 
of vortices and their size. The column density of CH + is there- 
fore highly dependent on «h (varies roughly as « H 2 ) while those 
of CH and W 2 {J > 3) are almost insensitive to the density in the 
range of parameters explored. Last, the turbulent rate of strain a 
has a weak effect: all the column densities decrease by less than 
a factor of 10 when a increases by a factor of 100. 

6.2. Submillimeter and millimeter observations: absorption 
lines towards continuum sources 

Conspicuous correlations have been observed in the diffuse ISM 
among the observed column densities of OH, H2O, HCO + , CO, 
C2H, CN, HCN and HNC seen in absorption along lines of sight 
towards extragalactic radio sources or star forming regions. The 
column densities of several species appear to be linearly corre- 
lated. Some are tight correlations such as N(OH)/N(HCO + ) = 
27 + 7 (Lucas & Liszt 1996), N(CN)/jV(HCN) = 6.8 + 1, and 
jV(HNC)/JV(HCN) = 0.21 + 0.05 (Liszt & Lucas 2001) in the 
direction of extragalactic radio sources. Others are looser corre- 
lations, still linear, such as N(U 2 0)/N(UCO + ) ~ 6 (Olofsson et 
al. submitted) and N(H 2 0)/N(OH) ~ 0.3 (Neufeld et al. 2002) 
found in several velocity components towards galactic star form- 
ing regions. The column densities of CO and C2H appear non- 
linearly correlated to those of HCO + (Liszt & Lucas 1998; Lucas 
& Liszt 2000). In all cases, these correlations are observed over 
a dynamic range of column densities as large as 30-100. 

To determine whether these large dynamic ranges corre- 
spond to actual variations of molecular relative abundances or 
are due to a large range of column densities of gas sampled, we 
have estimated the total amount of gas Nn = N(Hl) + 2N(H 2 ) 
along the lines of sight used in our study, whenever it was pos- 
sible. Using the A21 cm observations of HI (Du Puy et al. 1969; 
Radhakrishnan et al. 1972; Lazareff 1975; Dickey et al. 1978) 
on the one side, and the measurement of the A9 cm line of CH 
and the remarkable correlation between CH and H2 in the dif- 
fuse ISM (Liszt & Lucas 2002) on the other side, we estimated 
the total H column density for several lines of sight diplayed in 
Fig.[9land[T0l We found that most lines of sight sample about 1 
magnitude of gas, with the exception of the upper point of pan- 
els (a,b) that corresponds to 2.5 x 10 21 < N n < 3.6 x 10 21 crrT 2 
and the lowest for which A^h = 6 x 10 20 cm 2 . The next 2 points 
both verify A^h > 1.2 x 10 21 crrT 2 . The dynamic range of 30-100 
observed on the column densities in Fig. |9"1 and [TOl is therefore 
due in a small part to the amount of gas on the line of sight, 
and the actual dynamic range of observed relative abundances, 
N(X)/Nn, remains as large as about 10 along these diffuse lines 
of sight. 

All these data are displayed in Fig.[9land[T0l The left panels 
show the predictions of several PDR models while the right pan- 
els display those of the TDR models for comparison (normalized 
to 1 magnitude of diffuse gas sampled on the line of sight). The 
calculations were made assuming a recombination rate of HCO + 
of k = 2.4 x 10- 7 (r/300K)-°- 69 cm" 3 s" 1 (Ganguli 1988). This 
recombination rate being critical for some of the molecules pre- 



sented here, we discuss the impact of this choice in Appendix 

6.2.1. OH, H z O and HCO + 

For these molecules, the PDR models predict linear correlations 
between their abundances, but they fail to reproduce the ob- 
served ratios and the observed dynamic range of absolute abun- 
dances. 

The TDR models, without fine-tuning of the parameters, are 
consistent with the data over a broad range of gas densities 10 
cm 4 < n H < 200 cm 4 and rates of strain 10~ n s _1 < a < 10~ 10 
s _1 . Better agreement is obtained for small values of the rate of 
strain, i.e. a chemistry dominated by ion-neutral drift heating. 
Not only are the abundance ratios correctly reproduced by the 
TDR models, but the dynamic ranges are also reproduced even 
with a single value of ~s bs- Moreover these results are obtained 
without producing water abundances in excess of observed val- 
ues. The models at high densities are consistent with the up- 
per limit reported in HST-GHRS observations by Spaans et al. 
(1998). 

6.2.2. C 2 H and HCO + 

In the case of C 2 H, we find that the observed column densi- 
ties are reproduced by several models, something that the PDR 
models cannot even approach. The average abundance ratio 
jV(C 2 H)/jV(HCO + ) = 14.5 + 6.7 (Lucas & Liszt 2000) is re- 
produced for almost the same models that match the observed 
column densities of OH and H 2 0: 30 cm~ 3 < «h < 200 cm" 3 
and rates of strain 10~ n s _1 < a < 10~ 10 s _1 . 

6.2.3. CO and HCO + 

In the case of CO, the TDR models are closer to the observed 
column densities than the PDR model predictions, and again, 
high densities and low rates of strain are more favorable. In 
Appendix iBl we show, however, that if the energy dissipated in 
a vortex lifetime ry is reduced by a factor of 10 (i.e. the im- 
portance of the relaxation phase is increased by a factor of 10), 
the whole range of observed CO column densities is naturally 
reproduced for /in > 100 cirT 3 without significantly modifying 
the results regarding the other molecules. 

The dumpiness and large fluctuations of density along the 
line of sight - often required to find an agreement between ob- 
served and predicted column densities in UV-driven chemical 
models (e.g. Black & Dalgarno, 1977; Le Petit et al. 2006) - are 
not required here. 

6.2.4. CN, HCN and HNC 

The TDR models fail to reproduce the high column densities 
observed for CN, HCN and HNC, although the results are one 
to two orders of magnitude above the PDR models predictions. 
This discrepancy might be due to the nitrogen chemistry which 
involves neutral-neutral reactions whose rates are poorly known 
(Pineau des Forets et al. 1990, Boger & Sternberg 2005). 

6.3. The departure of carbon from ionization equilibrium 

Fig. QT| shows that the neutral carbon abundance predicted by 
the TDR models is higher by up to a factor of 10 than that of the 
PDR models where carbon is in ionization balance. This may 
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be related to the finding of Fitzpatrick & Spitzer (1997) that the 
electron density inferred from C and C + is higher than those in- 
ferred from other pairs of ions and neutrals (Mg, S, Ca). 




o 1 — 1 — 1 — ■ — 1 — ■ — ■ — ■ — ■ — 1 — ■ — ■ — ■ — ■ — 1 — ■ — ■ — ■ — ■ — 1 

50 100 150 200 

n H (cm" 3 ) 



Figure 11. Column densities ratio N{C)/N(C + ) computed with 
the PDR and TDR codes, as functions of the density n H . 



7. Discussion 

7.1. The gas homogeneity 

In our TDR models, the ambient chemistry is not treated with 
accurate UV radiative transfer. Instead, we assume that all the 
CNM on the line of sight has the same shielding from the am- 
bient UV radiation field, no matter its density and column den- 
sity. This takes into account the fractal structure of the CNM 
(Elmegreen et al. 2001) also found in the numerical simulations 
of Audit & Hennebelle (2005) and similar to that of edges of 
molecular clouds (Falgarone et al. 1991) and translucent molec- 
ular clouds (Stutzki et al. 1998). It also takes into account the fact 
that fragments that bear molecules are not isolated in space but 
are shielded from the ambient UV field by their neighbouring 
fragments (see Dobbs, Glover & Klessen 2008). This however 
does not have a major impact on our results for all the species 
that form predominantly in the warm chemistry, i.e. most of the 
species discussed here. 

7.2. The assumption of isochoric evolution 

During the vortex stage, the incompressibilty is reasonably jus- 
tified. (1) Radially, the advection force is 3 to 10 3 times larger 
than the neutral thermal pressure gradients V,. (n„kT„) over the 
whole structure. Since it is also 10 2 to 10 3 times greater than the 
friction force Fj„ (see Sect. 12.31) . V,- (n„kT n ) ~ 1 - 30F,„. The ion 
thermal pressure gradients V,. (rtjkTi) ~ 2 ■ 10~ 4 V,. (n n kT n ) are 
consequently negligible compared to F The assumption in Eq. 
( fTTb . involved in the computation of the steady state, is therefore 
jusitified. (2) Along the z axis V z (n„kT„) is smaller than the ad- 
vection force as long as a > 10~ 10 s _1 . For a < 10~ 10 s these 
terms are comparable. In this case, the axial expansion timescale 
is ~ 1000 yr, within a factor of 3 of the vortex lifetime (see Fig. 

During the relaxation phase, the thermal pressure gradi- 
ents are dominant and therefore drive its dynamical evolution. 



According to the temperature gradients of Fig. [2] (panel c), the 
inner regions (r < 0.5ro) should be compressed while the outer 
regions (r > 1.5ro) should expand. The characteristic timescales 
associated with those motions are found to be short compared to 
the chemical timescales (see Sect. 14.3b . However, an isochoric 
relaxation was assumed because the regions that are chemically 
the richest and have the largest contribution to the column densi- 
ties are those where the thermal pressure gradients are the small- 
est (0.5r < r < 1.5r ). 



7.3. Loosening the vortex framework 

The vortex is a model that has the great and unique advantage of 
analytically coupling the large and small scales, but the knowl- 
edge we have of its lifetime is drawn from laboratory experi- 
ments only. This lifetime Ty however plays an important role in 
the modelling of a diffuse line of sight: Eq. d23l shows that if a 
species is predominantly produced during the relaxation phase, 
its column density is proportional to 1 fry. 

We have made the choice in Sect. 15.2.31 to infer t v from the 
total energy E dissipated in the vortex, under the constraint that 
Ty remains larger than the vortex period. We show in Appendix 
|B]how the TDR model results are modified when E (and there- 
fore Ty) is ten times larger or smaller than the adopted value. As 
expected, the impact of E is greater for the molecules that have 
a long relaxation time. Fig. IB. 11 shows that a small value of E 
leads to a better agreement with the available observations. This 
is shown for OH and CO but the same behavior is observed for 
CN, HCN and HNC. Such a result is in favour of very short 
bursts of turbulent dissipation, for which Ty may become (at 
small a) shorter than the vortex period. 

A non steady state model of a short-lived but intense 
velocity-shear may thus be more realistic and produce an equiv- 
alent chemical enrichment. This suggests that the vortex frame- 
work is not essential. 



7.4. The chemistry 

One source of uncertainties is undoubtedly linked to the chem- 
istry itself. While the results are weakly dependent on the dy- 
namic parameters a and rtu, a few reaction rates are critical. This 
is illustrated by Fig. lA.ll that displays the predictions of the PDR 
and TDR models for a recent value (measured in the laboratory, 
Mitchell & Mitchell 2006) of the recombination rate of HCO + . 

The uncertainties on other chemical rates such as the recom- 
bination rate of HCNH+ (Mitchell & Mitchell 2006) and the 
CN photodissociation rate (Kopp 1996), and neglected reactions 
(with negative ions for instance, see Dalgarno & MacCray 1973) 
might affect our modelling. 

In this respect, a promising route would be to further link 
the observations of Sect. 16. H and those of Sect. |6.2| by using the 
species that are detectable in both UV and radio domains: for 
instance CH, CN, OH and CO. These molecules would permit us 
to establish valuable correlations between CH + and the oxygen 
bearing species. It would also be a way to evaluate the column 
density of molecular hydrogen using two independent methods: 
the remarkable correlation between CH and H2 (Liszt & Lucas 
2002) and the less reliable correlation between CN and H2 (Liszt 
& Lucas 2001). 
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8. Conclusions and perspectives 

We have built TDR models of diffuse gas in which the gas tem- 
perature and chemistry are driven by small-scale bursts of turbu- 
lence dissipation. We use the framework of a modified Burgers 
vortex to analytically couple the large scales of the ambient tur- 
bulence to the small scales where dissipation actually occurs, 
and to compute the ion-neutral drift generated by the large neu- 
tral accelerations in the vortex. The resulting timescales are short 
and comparable to those of chemical evolution, which necessi- 
tates a non-equilibrium thermal and chemical approach. 

The main feature of these TDR models is that, for the first 
time, we quantify the coexistence on a random line of sight 
across the medium, of a number of vortices in a stage of active 
dissipation with gas in thermal and chemical relaxation, after 
the end of the dissipation burst. We also constrain the number of 
dissipative structures on a line of sight by the average turbulent 
energy available in the local ISM and its transfer rate in the cas- 
cade. The key parameter is the turbulent rate of strain a due to 
the ambient turbulence. 

We find that these bursts of dissipation, short-lived and local- 
ized, fill at most a few percent of a random line of sight but have 
a measurable impact on the molecular abundances in the diffuse 
medium. For a broad range of rates of strain and densities, the 
TDR models reproduce the CH + column densities observed in 
the diffuse medium and their correlation with highly excited H2. 
They do so without producing an excess of CH. 

As a natural consequence, they reproduce the correct abun- 
dance ratios of HCO+/OH and HCO + /H 2 0, and the dynamic 
range of about one order of magnitude over which they are ob- 
served. Larger C2H and CO abundances than found in other 
types of models, are additional outcomes of the TDR models 
that compare reasonably well with the observed values and their 
relation to the HCO + abundances. Those results are found for a 
broad range of physical parameters, with rates of strain in the 
range 1CT 11 s" 1 < a < 1(T 10 s _1 . 

We find that neutral carbon exceeds the abundance expected 
at ionization equilibrium, in agreement with fine-structure line 
observations. The abundances and column densities computed 
for CN, HCN and HNC are one order of magnitude above PDR 
predictions and close to the observed values, although a discrep- 
ancy still exists. 

The comparison with observed column densities favors 
chemical enrichment dominated by ion-neutral friction, involv- 
ing shear structures of radius ~ 100 AU or more. It also favors 
short dissipation bursts. Most of the species are then produced 
during the relaxation phase. In this case, their line profile loses 
the dynamic signatures of the vortex. The fact that some species 
are formed in the active phases and some others in the relaxation 
phases might explain the discrepancy among the physical gas 
parameters inferred from different species along similar lines of 
sight. 

It is foreseeable that the vortex framework and the assump- 
tion of an isochoric relaxation will have to be superseded by nu- 
merical simulations to: 

(1) compute the decoupling of magnetic field and neutrals in the 
intense small-scale velocity-shears generated by intermittency, 
with boundary conditions imposed by large scale turbulence, 

(2) take into account the effect of the field on the velocity shears, 

(3) accurately involve all the pressure gradients (thermal and 
magnetic) in the relaxation phase, and 

(4) take advantge of the huge amount of information contained 
in the shape of the line profiles. The difficulty will remain to 
couple the large and small scales, over at least 5 orders of mag- 



nitude, that is critical to satisfy the energy requirements of the 
TDR models and is currently beyond the numerical capabilities. 
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Later the same group (Mitchell & Mitchell 2006) revised their previous value by 
taking into account the exploration of the low energy domain (< 0.01 eV). They 
found : 

k = 0.7 X 10" 7 X (7//300KT 5 cmV . (A.3) 

More recently Korolov et al. (in press, flowing aftergow plasma) performed mea- 
surements in the range 150K < T < 270K. They obtained results in good agree- 
ment with those of Le Padellec et al. (1997), this is to say a steep dependence on 
the temperature: 

k = 2.0 X 10" 7 X (7//3O0KT 1 3 em's" 1 . (A.4) 

Fig. lA. ll displavs the results we obtain with the TDR and PDR models as- 
suming the rate of Mitchell & Mitchell (2006) which is the smallest yet found. In 
comparison Fig.|9]was obtained assuming the rate of Ganguli et al. (1988) which 
is one of the highest and is usually adopted in other chemical networks (UMIST 
database; OSU database). 

Appendix B: Influence of the total energy E 
dissipated per vortex. 

Fig. |B.l| illustrates the effect on the column densities of OH, CO and HCO + of 
varying the energy dissipated in the vortex E by two orders of magnitude. 

Fig. IB. II suggests that small values of E lead to better agreement with the 
observations. Smaller amounts of dissipated energy, in our framework, means 
shorter vortex lifetimes, that accordingly give a larger relative importance to the 
relaxation phase (Eq. |23) . In the models displayed on the top panels, the life- 
time is such that the relaxation phase dominates the production of almost all the 
species we are interested in, namely CH, C 2 H, OH, H 2 0, HCO+, CO, CN, HCN 
and HNC. 

Appendix C: UV-dominated versus 

turbulence-dominated chemical networks 

Figs. IC.ll and |C.2| display the main production and destruction routes of the 
molecules of interest in the ambient UV-dominated diffuse medium and in the 
vortex for the reference TDR model at a radius r = tq respectively. 

These figures are simplified: only the dominant reaction of production and 
the dominant reaction of destruction of the species are shown. For example one 
would expect the reaction O + Hj — > OH + + H 2 to be displayed in Fig. IC.ll 
But in the model presented here, the efficiency of this reaction on the production 
of OH+ is only 12%, while the efficiency of + + H 2 -> OH+ + H is 86%. 
Therefore it does not appeal' on the diagram. 
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Appendix A: The recombination of HCO + 

One critical reaction rate of our work is the recombination of HCO + which de- 
termines the equilibrium between the oxygen bearing molecules (OH, CO, H 2 0, 
...). The associated rate k has varied over more than an order of magnitude during 
the last 30 years among various experiments. 

Adams et al. (1982) and Amano (1990) performed measurements of k 
(studying respectively a flowing and a stationary afterglow plasma) in a tem- 
perature range of 100K < T < 300K and found respectively 1.1 X 10" 7 cm 3 s _1 
at 300 K and 3.1 X 10" 7 cm 3 s" 1 at 273 K. In 1988 Ganguli et al. (stationary 
afterglow plasma) explored the high energy domain 293K < T < 5500K and 
obtained a scaling law: 

k = 2.4 X 10" 7 X (77300K)" - 69 cmV 1 . (A.l) 



In approximatively the same range of energy (but with the merged beam mea- 
surement technique) Le Padellec et al. (1997) obtained 



k = 1.7 X 10" 7 (T/300K)"' 2 em's" 1 . 



(A.2) 
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Figure 8. Observations compared to PDR models (left panels) and to TDR models (right panels). Data (open circles) - CH column 
densities are from Crane et al. (1995) and Gredel (1997). CH + column densities are from Crane et al. (1995), Gredel (1997) and 
Weselak et al. (2008). H* 2J=4 column densities are from Spitzer et al. (1974), Snow (1976, 1977), Frisch (1980), Frisch & Jura (1980) 
and Lambert & Danks (1986). The excitation diagrams of H2 are from Gry et al. (2002) and Lacour et al. (2005). The data on panels 
(e,f) and (g,h) are scaled to A^h = 1.8 x 10 21 crrT 2 and A^(H2) = 5 10 20 cirT 2 respectively, and correspond to lines of sight toward 
stars of type later than B2. For clarity, on panels (g,h), the points for a given level are slightly shifted on the abscissa. PDR and 
TDR models (filled symbols) - Computed for several densities: 10 (squares), 30 (triangles), 50 (crosses) and 100 cirT 3 (circles). 
All the TDR models are computed for Ay = 0.2. The models of panels (e,f) and (g,h) have been scaled as the data. In panel (f) the 
TDR models are computed for a varying along each curve between 10" (top right) and 10~ 9 s~' (bottom left). 
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Figure A.l. Observations compared to PDR models (left panels) and to TDR models (right panels). Data (open circles) - Same 
as Fig. [9] PDR and TDR models (filled symbols) - Same as Fig. [9] except for the recombination rate of HCO + chosen to be 
0.7 x 1(T 7 (77300 KT - 5 cm 3 s" 1 (Mitchell & Mitchell 2006). 



I 

E 
X 

O 



B. Godard et al.: Models of turbulent dissipation regions 

E/10 E/10 

1e+17 w 



21 



1e+15 f 



1e+14 r 



1e+13 r 



"1 0-1-1 2 i ■ i i i i 

1e+10 1e+11 1e+12 1e+13 

-2x 




E 
O 

o 



- 1 — I I I I I III 1 — I I I I Mil 1 — I I I I 111 




N(HCO + ) (cm"') 



1e+16 r 
1e+15 r 
1e+14 r 



"|0-|-13 * 1 I i I i 

1e+10 1e+11 1e+12 1e+13 

-2x 



N(HCO + ) (cm"') 



I 

E 

X 

O 



1e+15 r 
1e+14 r 
1e+13 r 
1e+12 - 



-i — i i i ■ i ■ ■ I 1 — i i 1 1 mi 1 — I I I I III 




J I I I I 



C\l 

E 
O 

o 



1e+10 1e+11 1e+12 1e+13 
N(HCO + ) (cm" 2 ) 



— i — i i ■ 1 1 ■ ■ | 1 — i i 1 1 1 ■ ■ | 1 — i i 1 1 in 

HSH H 

t 



Wffc leaf 



1e+17 
1e+16 
1e+15 
1e+14 L 



"|0-|-13 I 1 I i I i 

1e+10 1e+11 1e+12 1e+13 

-2x 





N(HCO + ) (cm"') 



Ex 10 



Ex 10 



<M 

E 
X 

O 



1e+15 r 
1e+14 r 
1e+13 r 
1e+12 - 



— I — I I I I INI 1 — I I I I Mil 1 — I I I I III 




j * i ' i 



CM 

E 

6 
o 



1e+10 1e+11 1e+12 1e+13 
N(HCO + ) (cm" 2 ) 



1e+17 r 
1e+16 r 
1e+15 r 
1e+14 r 
1e+13 '- 



— i — i i 1 1 1 ■ ■ | 1 — i i 1 1 1 ■ ■ | 1 — i i 1 1 in 




_i i i i i 



1e+10 1e+11 1e+12 1e+13 
N(HCO + ) (cm" 2 ) 



Figure B.l. Observations compared to the TDR models. The data are from Lucas & Liszt (1996) and Liszt & Lucas (1998). The 
TDR models are scaled to A^h = 1.8 x 10 21 cirT 2 and are computed for several values of the total energy E dissipated by a vortex: 
from 0.1 (top panels) to 10 times (bottom panels) that of the reference model. The symbols are the same as in Fig. [9] 
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1 Carbon Carbon Carbon & O^gen O^gen 

Figure C.l. Chemical network of a UV-dominated chemistry: «h = 30 cm -3 and Ay = 0.4. This figure is simplified: for each species 
only the dominant reaction of production and the dominant reaction of destruction are displayed. 




2. Carbon Carbon Carbon & O^gen 0*>{gen 

Figure C.2. Same as Fig. lC.ll for a turbulence-dominated chemistry: «h = 30 cirT 3 , A y = 0.4, a = 3 x 10 _1 1 s _1 at a radius r — ro. 
The red arrows show the endoenergetic reactions with the energy involved. 



